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ABSTRACT 


A  comparative  study  between  GASBOGOY  and  (HOME  data  shows 
considerable  differences  in  observed  travel  times,  amplitudes, 
and  in  the  excitation  of  surface  waves.  Die  Pn  and  P  arrivals 
from  GASBGGOY  follow  the  Jeffreys-Bullen  curve  closely,  whereas 
the  arrival  times  for  GNOME  were  early  to  the  east  and  late  to 
the  west.  Relatively  higher  amplitudes  of  phases  associated 
"ith  Fig  and  Lg  were  recorded  to  the  east  of  GASHJQGY,  whereas 
Pn  and  P  and  surface  waves  (LR)  recorded  higher  amplitudes  to 
the  west  of  OASBOGGY.  Rayleigh  waves  were  strongly  recorded 
at  large  distances  from  QASBQQGY  with  an  average  magnitude  of 
U.09,  as  compared  with  a  Mg  and  Mg  of  £>.21  and  U;73»  respectively. 
Since  P  and  PKP  data  were  not  recorded  beyond  86  degrees,  the 
body  wave  magnitude  computed  according  to  Evernden  (Mg)  appeared 
to  be  a  more  reasonable  value;  A  refinement  of  the  QASBDQOY 
location  was  attempted  by  adjusting  the  arrival  times  relative 
to  the  residuals  observed  from  GNCME.  The  maximum  epicentre! 
error  increased  approximately  7  kilometers  to  the  north  and  22 
kilonaters  to.  the  west. 
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INTRODUCTION 


ftroject  GASBUQQY,  the  first  cooperative  Government-industry 
Plowshare  project,  was  detonated  10  December  196?  at  19h  30m 
00. le  U.T.  at  a  depth  of  1*21*0  feet  in  a  shale  environment.  The 
26-kiloton  nuclear  explosion  was  located  in  the  San  Juan  Basin 
near  Farming ton,  New  Mexico.  The  geographio  coordinates  were 
36°l*0'l*0.i*"  N,  107°i2'30.3"  W.  The  primary  purpose  of  the 
experiment  was  to  test  a  proposed  method  of  increasing  the  pro¬ 
duction  of  oil  and  gas. 

In  this  report  the  results  of  GASBUGGY  will  be  compared 
with  data  from  GNOME,  an  underground  explosion  detonated  in  a 
salt  mine  near  Carlsbad,  New  Mexico.  Apparently  due  to  its 
location  on  the  boundary  between  the  Great  Plains  and  the  Rocky 
Mountains,  GNGME  exhibited  considerable  differences  between  ampli¬ 
tude  attenuation  and  seismic  velocities  associated  with  wave 
propagation  to  the  east  and  to  the  west  of  the  test  site.  In 
contrast^  the  QASBUOGY  site,  which  was  located  approximately 
l*°N  and  3°W  of  GNOME,  lies  well  within  the  Rocky  Mountain  region. 
Che  of  the  principal  objectives  of  this  report  will  be  to  deter¬ 
mine  the  effect  of  shifting  the  test  site  several  hundred  kilo¬ 
meters  to  the  west  in  a  region  where  rapid  changes  in  crustal 
and  upper  mantle  structure  are  suspected. 
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The  GASBGQOY  test  site  was  located  in  the  east-central 
portion noY  the  San  Juan  Basin.  The  basin,  approximately  90 
miles  in  diameter,  an  area  of  about  10,600  adits,  has  an  esti¬ 
mated  6000  feet  of  structural  relief.  Otologic  formations 
exceeding  15,000  feet  in  total  thickness  renge  from  pre-Cambrian 
age  to  Eocene  in  the  interior  part  of  the  basin.  The  basin  is 
rimmed  by  monoclines  on  the  north,  east,  and  vest.  The  rocks  in 
the  southern  part  of  the  basin  indicate  periods  of  diastrophism, 
including  folding,  faulting,  mid  igneous  activity. 

DISCP5SICH  OF  DATA 

GASBQQGY  travel  times  were  available  from  157  stations,  and 
ground  amplitudes  were  determined  at  129  stations  where  instru¬ 
mental  responses  are  veil  known.  The  closest  station  was  at  a 
distance  of  liu  kilometers  from  the  test  site,  and  the  most 
distant  was  9U67  kilometers  at  Matsu  a  hiro,  Japan.  Seismic 
data  were  recorded  from  19  stations  outside  of  the  North 
American  Continent- -two  each  in  Greenland,  South  America, 

Finland,  SVeden,  and  Japan*  and  nine  in  Europe. 

Much  of  the  travel -tims  information*  crams  from  the  World-Wide 
stand*—’  seismograph  stations  and  other  similar  permanent  obser¬ 
vatories  operated  by  universities  and  the  Coast  and  Qeodetic 

Survey.  At  distances  less  than  500  kilometers,  ll*  temporary 
stations  operated  by  the  Ooast  and  Geodetic  Survey  and  the 


Geological  Survey  provided  important  supplemental  travel  times. 

In  this  report,  only  ten  stations  contributed  useable  data.  The 
temporary  station  operational  data  are  summarized  in  Appendix  I. 

The  temporary  stations  were  equipped  with  high-  and  low-gain^ 
three -component  ASCI*  and  ASC2*  Seismic  Data  Systems.  To  aid  in 
the  comparison  with  the  GNOME  event,  LRSM  data  are  included  in 
the  various  diagrams.  Ihe  LRSM  stations  were  equipped  with  three- 
component  short-  and  long-period  instruments.  The  horizontals 
were  oriented  radially  and  transverse  to  the  great  circle  path 
from  the  QASHJQGY  epicenter. 

The  basic  data  used  In  this  report  are  summarized  In  Figures 
1,  2,  and  In  Table  1.  Figure  1  shows  the  geographic  location  of 
North  American  stations  and  type  of  signal  reception  at  the  stations 
known  to  have  recorded  QASHJQGY.  Figure  2  shows  reception  character 
istics  at  other  teleseismlc  stations.  In' Table  1,  the  basic  data 
Includes  distance  (degrees  and  kilometers),  Instrumental  component, 
phase  identification,  arrival  time  (uncorrected  for  ellipticity 
and  station  elevation),  period  and  representative  amplitude  of 
selected  phases,  and  both  Oaten  berg -Richter  body-wave  magnitude 
and  the  magnitude  defined  by  Evernden  (1967)  which  relates 
regional  apparent  velocities  with  representative  amplitude 
attenuation.  The  Jeffreys- Bullen  residual  (observed-computed) 
takes  into  account  the  ellipticity  and  elevation  correction  for 
the  longitudinal  phases.  Figure  3  illustrates  signal  detection 

*Fortable  magnetic  tape  systems  designed  and  built  by  Albuquerque 
Seismological  Center  (C&GS) 
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as  a  function  of  instrumental  magnification  and  distance.  As 
seen  previously,  there  was  considerable  overlap  of  reception 
and  nonreception  at  similar  distances,  which  indicates  signal 
detection  is  not  a  simple  function  of  range  and  instrumental 
magnification.  Other  factors,  such  as -the  signal -to -noise  ratio 
and  the  pattern  of  regional  amplitude  anomalies,  appear  to  play 
an  important  role. 

In  the  accompanying  graphs  of  the  data,  which  correspond 
to  the  phases  fti  and  P,  Pg,  Ig,  and  LR,  the  stations  to  the 
east  of  the  test  site  are  separated  from  those  to  the  west  by 
different  symbols.  This  separation  roughly  approximates  the 
division  of  the  Continental  United  States  into  two  major 
seismic  provinces  as  delineated  by  the  Rocky  Mountain  Rront. 

The  R“3  reference  attenuation  line  was  visually  fitted  to 
the  data. 

AMPLITUDES  -  TRAVEL  TIMES 
Ph  and  P 

Amplitudes  of  the  largest  motion  in  the  first  three  cycles 
are  illustrated  in  Figure  h.  The  attenuation  rate  of  the  waves 
appeared  to  fit  the  R~^  reference  line  out  to  7OU-0OO  kilometers 
as  has  been  observed  in  the  case  of  the  NTS  events.  Beyond  this 
point  there  was  no  apparent  correlation,  of  course,  as  phases 
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representing  uaeper  penetration  became  first  arrivals. 

As  will  be  seen  in  this  report,  estimates  of  amplitude 
attenuation  rates  at  regional  distances  consist  of  (l)  matching 
the  identified  phase  with  an  appropriate  apparent  velocity;  (2) 
determining  the  amplitude  decay  associated  with  the  corresponding 
refractor. 

Several  differences  were  noted  when  GNOME  and  GASHJGGY 
data  were  compared.  The  higher  amplitudes  of  the  eastern 
stations  at  similar  distances  to  those  to  the  west  as  observed 
for  GNGKE  were  not  reflected  by  the  GASBUGGY  data.  Generally, 
the  stations  to  the  west  tended  to  have  equal  or  higher  ampli¬ 
tudes. 

As  shown  in  Figure  5»  the  arrival  times  of  GASBUGGY  followed 
the  Jeffreys -Bull en  curve  ciosely.  GNCME,  in  the  distance  range 
of  1000  kilometers,  reported  eastern  stations  about  5  sec  early, 
as  compared  with  3  sec  for  GASBUGGY.  Arrival  times  to  the  west 
of  GASBUGGY  showed  a  more  striking  disagreement  in  the  distance 
range  of  1000  to  1$00  kilometers.  GNOME  reported  a  maxi  mum 
residual  of  +1  sec  whereas  GASBUGGY  reported  a  maximum  of  +2.5 
sec.  Also,  a  significant  number  of  stations  to  the  west  of 
GASBUGGY  reported  early  arrival  times,  whereas  stations  west 
of  GNOME  generally  reported  late  arrivals. 

At  distances  beyond  20  degrees,  the  average  recorded  arrival 
times  were  about  1.0  sec  earlier  than  the  travel  times  predicted 
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by  the  Jeffreys -Bullen  table.  Early  arrival  times  were  also 
observed  from  the  GNOME  explosion,  as  well  as  from  earlier  ex¬ 
plosions  in  the  Pacific  area  (Gutenberg  and  Richter,  196b, 

Burke -Gaffney  and  Bullen,  1957;  Carder  and  Bailey,  1958;  Carder, 
Gordon,  and  Jordan,  1966). 

Travel -time  data  were  available  in  the  distance  range  from 
Ui  to  9500  kilometers.  All  stations  closer  than  1*00  kilometers 
were  located  along  a  rougi  profile  south  of  the  test  site.  A 
straight-line  approximation  of  the  first  leg  (t^)  of  the  t  ravel - 
time  curve,  corresponding  to  the  arrival  at  station  UH  at  a 
range  of  U h  kilometers,  yields  an  apparent  velocity  of  5 - 3U 
km/ sec.  The  second  branch  of  the  travel -time  curve,  with  an 
apparent  velocity  of  6.29  km/sec  and  a  time  axis  intercept  of 
l.U  sec,  corresponds  to  the  observed  first  arrivals  covering 
the  50  to  circa  180  kilometers.  Several  investigations,  in¬ 
cluding  Carder,  et.  al. ,  (1966)  and  Evemden  (1967),  have  indi¬ 
cated  that  sub-Moho  refractors,  representing  apparent  velocities 
of  7.9,  8.5,  and  10.5  km/sec,  exist  beneath  the  western  United 
States.  Reference  lines  representing  these  apparent  velocities 
were  visually  fitted  to  the  station  data  shown  in  Figure  5.  As 
can  be  seen  in  Figure  5,  there  was  considerable  latitude,  espe¬ 
cially  at  the  crossover  points,  in  placing  the  reference  lines 
through  the  data  points.  Possibly,  the  various  apparent  velocities 
associated  with  the  stations  recording  first  arrivals  were 
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dependent  upon  the  signal- to-noise  amplitude  ratio.  In  addition 
to  the  graphic  representation  in  Figure  5  all  travel  times  are 
presented  in  Table  2. 

TABLE  2  —  Summary  of  Travel  Times 


Surface  Crossover 


Travel  Time 


Intercept 


tx  -  A/5.3U 

t2  -  A^6.29 

+  l.U 

50 

t3  -  A/7.9 

♦  7.0 

155 

tU  -  4^8.5 

+  15.8 

910 

Regionalization  was  apparent  when  the  station  travel-time 
residuals  were  plotted  and  contoured  on  a  map  (Figure  6).  Hie 
region  of  greatest  negative  residuals  was  located  in  the 
Mississippi  River  Valley,  and  the  area  of  greatest  positive 
residuals  was  located  in  the  California  area.  This  latter 
feature,  although  not  as  well  defined  as  with  GNCME,  could 
represent  the  effect  of  the  Sierra  Nevada  "root”  projected  on 
the  surface  to  the  west.  The  map  generally  agreed  with  the 
GNCME  map,  with  the  exception  that  GASBUOGY  indicated  negative 
residuals  located  in  the  Utah-Nevada  area. 

Figure  7  illustrates  the  patterns  of  maximum  amplitudes 
of  Pn  and  P  waves  across  the  United  ^ates.  The  contour  interval 
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was  based  on  the  geometric  progression  1,  2,  k,  8,  16,  32,  etc., 
to  permit  contouring  over  the  large  range  of  observed  amplitudes. 
Although  the  distance  range  was  different  for  the  two  events, 
several  geographic  regions  showed  similar  amplitude  anomalies 
for  GNOME  and  GASBUGGY.  For  example,  the  low-amplitude  anomaly 
in  Arkansas,  at  a  distance  of  1000  kilometers  from  GNOME,  was 
repeated  for  GASBUGGY  even  though  the  distance  was  1200  to  11*00 
kilometers.  The  high  anomaly  which  occurred  in  Texas,  600  kilo¬ 
meters  from  GNOME,  was  not  as  well  delineated  by  GASBUGGY.  However, 
there  was  less  control  for  GASBUGGY.  The  low -amplitude  area  which 
occurred  in  northern  Colorado  from  GNOME  had  shifted  northward 
from  GASBUGGY  but  had  still  maintained  the  1000 -kilometer  spacing. 
Also,  as  found  in  GNOME,  GASBUGGY  indicated  low-amplitude  anomalies 
in  the  California  area  and  around  the  Nevada  Test  Site.  TVo  high- 
amplitude  anomalies  were  shown,  due  to  greater  station  control  by 
the  recorded  GASBUGGY  data:  One  was  2800  kilometers  to  the 
northeast  and  the  other  1700  kilometers  to  the  northwest. 

In  Figure  7,  the  apparent  velocities  of  the  station  first 
arrivals  are  indicated  above  the  station  code.  At  1000  kilo¬ 
meters  to  the  east  of  the  GASBUGGY  test  site,  the  ^7.9  velocity 
arrival  was  replaced  by  the  arrival  of  the  P0  ^  refractor,  and 
beyond  2000  kilometers,  only  the  arrival  with  an  app«w?nt  ve¬ 
locity  of  P  £  was  recorded.  To  the  west,  the  distinction 


between  the  P„  „  and  P  apparent  velocity  arrivals  were  not  as 
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abrupt  as  it  was  to  the  east.  In  the  distance  range  of  500  to 
1000  kilometers,  the  identification  of  the  Py#^  velocity  refractor 
was  probably  dependent  upon  a  favorable  signal-to-noise  amplitude 
ratio  at  the  individual  stations.  The  appearance  of  the  Py^ 
arrival,  as  suggested  by  Evemden  (196?),  could  be  distinctive 
in  differentiating  between  the  Eastern  and  Western  United  States. 
With  the  exception  of  WMO  and  GV-,  the  Rocky  Mountain  Front  was 
roughly  parallel  to  the  alinement  of  stations  which  recorded  the 
Py  ^  velocity  arrival.  As  contrasted  to  the  east  where  the 
P10.  £  velocity  arrivals  were  observed  only  beyond  2000  kilometers, 
the  P10.5  arrival  first  appeared  at  a  range  of  1800  kilometers 
to  the  northwest. 

lhe  preliminary  conclusion  based  on  the  station  apparent 
velocity  arrivals  shows  that  the  complexity  of  the  upper  mantle 
structure  results  in  significant  changes  in  the  P-phase  arrival? 
that  cannot  be  accounted  for  by  distance  alone.  As  illustrated 
in  Figure  7,  the  mantle  structure  is  more  complex  toward  the 
west  as  compared  with  the  Eastern  United  States. 

It  was  of  interest  to  note  that  at  ROD  and  WM-,  the  same 
distance,  velocity  refractor,  and  quadrant  had  striking  differ¬ 
ences  in  amplitudes.  The  low  amplitude  recorded  at  RCD  could 
possibly  be  the  result  of  a  shadow  effect  associated  with  the 
Black  Hills  physiographic  province. 


p& 

As  illustrated  in  Figure  8,  Pg  was  evident  out  to  a  distance 
of  I63O  kilometers  to  the  east  and  lliOO  kilometers  to  the  west. 
The  Pg  amplitudes  were  generally  higher  toward  the  east  than  to 
the  west,  and  the  wave  attenuation  could  be  represented  by  the  R“^ 
reference  line  out  to  1000  kilometers,  followed  by  a  more  rapid 
decay. 

l£ 

As  shown  by  Figure  9,  the  wave  attenuation  appeared  to 

follow  the  R~3  reference  line.  Die  Ig  waves  were  recorded  3975 

kilometers  to  the  east  and  only  2911  kilometers  to  the  west. 

Unlike  GNOME,  the  stations  to  the  east  had  higher  anplitudes  than 

those  to  the  west.  As  with  GNOME,  the  attenuation  rate  was  less 

rapid  toward  the  east  than  to  the  west. 

It  was  of  interest  to  note  that  Ig  waves  traversed  the 

Hudson  Bay  and  other  shallow  waters  south  of  RES.  This  was 

indicated  by  the  recorded  Ig  waves  at  both  RES  and  FBC.  This 

observation  corresponded  to  the  findings  of  Herrin  (i960). 

Although  a  strong  P  signal  was  observed  at  MBC,  Ig  waves  were 

« 

not  recorded.  This  verified  the  presence  of  an  oceanic  segment 
between  the  mainland  and  MBC. 

Unlike  GNOME,  Love  waves  were  recorded  on  WWNSS  and  LRSM 
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long-period  instruments  at  maximum  ranges  of  2665  kilometers  to 
the  east  and  1363  kilometers  to  the  west.  The  period  of  Love 
waves  ranged  from  8  to  30  sec,  with  the  dominant  period  about 
15  sec.  The  horizontal  component  which  recorded  the  wave  form 
best  is  indicated  by  a  special  symbol  in  Figure  10.  The  wave 
attenuation  rate  appeared  to  be  more  rapid  to  the  east  than  to 
the  west. 

LR 

As  observed  from  Figure  11,  Rayleigh  waves  generated  by 
GASBUGGY  were  recorded  beyond  UOOO  kilometers  to  the  east  and 
2000  kilometers  to  the  west,  whereas  GNCME  reported  Rayleigh 
waves  out  to  only  500  kilometers.  The  wave  form  designated  LR 
had  a  duration  of  several  cycles  with  a  period  of  approximately 
12  sec.  At  similar  distances,  in  the  distance  range  of  900  to 
2000  kilometers,  the  stations  to  the  west  had  higher  amplitudes 
than  the  stations  to  the  east. 

BODY  WAVE  MAGNITUDES 

Body  wave  magnitudes,  corresponding  to  the  P  wave  recorded 
by  short-period  vertical  instruments,  were  computed  according  to 
the  method  given  by  Gutenberg  and  Richter  (1956)  and  are  shown 
in  Figure  12  <  Note  the  scattering  of  data  when  the  amplitude - 
distance  relationships  were  applied  to  the  calculation  of 
body  wave  magnitudes.  Hie  average  magnitude  representing 
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stations  at  ranges  less  than  16  degrees  is  5.20  and  greater  than 
16  degrees  the  magnitude  is  5*03.  The  value  of  5*21  represents 
the  overall  average  magnitudes. 

A  second  plot  of  body  wave  magnitudes  (Figure  13)  has  been 
prepared  to  show  the  effects  of  regional  corrections  suggested 
by  Evemden  (1967).  In  this  plot,  magnitudes  in  the  range  of 
200  to  3000  kilometers  were  recomputed  using  "Q"  factors  which 
were  generated  by  the  formula  Mg  =  Q  +  q  +  S,  with  reference  to 
an  event  of  magnitude  5,  are  listed  in  Appendix  2. 

An  average  of  all  data,  covering  the  entire  distance 
range,  gives  a  magnitude  of  U .73 .  Stations  beyond  16  degrees 
indicate  a  value  of  In 99,  whereas  In 52  represents  the  data 
obtained  at  less  than  16  degrees.  The  application  of  the  new 
set  of  "Q"  factors  did  not  completely  remove  the  scatter  in 
the  magnitude  data.  In  some  instances,  especially  in  the  range 
of  1700  to  3000  kilometers,  the  adjusted  reductions  changed  the 
magni+'ide  residuals,  either  raising  or  lowering  by  several 
tenths,  depending  on  the  refractor  velocity.  Evidently  there 
are  other  factors,  such  as  radiation  patterns  and  ground  factors 
at  the  receiver,  which  contribute  significantly  to  the  data 
scatter.  However,  if  we  were  to  choose  between  the  first 
computed  value  of  5.21  and  the  adjusted  7alue  of  U.73»  we  would 
favor  U - 73  since  P  data  were  not  observed  beyond  86  degrees. 
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Generally,  an  "average”  5*2  event  would  generate  diffracted  P 
data  which  could  easily  be  read  on  WWNSS  seismograms. 

SURFACE  WAVE  MAGNITUDES 

The  Wbrld-wide  standard  stations  and  the  LRSM  stations, 
which  operate  three-component  long -period  systjms,  provided  the 
major  source  of  the  surface  wave  information.  Data  from  other  seis¬ 
mograph  stations  with  similar  instrumentation  were  used  whenever 
possible. 

Surface  wave  magnitudes  were  computed  with  the  use  of  a  formula 

Ms  -  log10A/T  +  1.66log10  A  °  +  3*3 

where: 

A  *  ground  amplitude  in  microns 
T  *  ground  period  in  seconds 
A-  distance  in  degrees 

developed  by  Vanek  (1962)  and  adopted  by  the  IASPEI  Committee 
which  recommended  that  the  formula  be  applied  to  20-sec  waves. 

This  was  not  possible  in  the  f;ase  of  GASHJGGY  because  the  only 
observable  waves  had  periods  of  circa  12  sec .  A/T  values  measured 
from  these  waves  were  entered  into  the  formula  to  obtain  Mg. 

Hie  Russians  have  indicated  that  the  formula  applies  to  other 
wave  periods  as  well  when,  the  measurements  of  A/T  are  associated 
with  maximum  amplitudes  (Soloviev  and  Shebalin,  1957;  Gutenberg, 
19U5). 

Based  on  the  formula,  the  computed  results  from  the  three - 
component  long-period  seismographs  are  presented  in  Table  3. 
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The  LRSM  stations,  as  contrasted  to  the  World-wide  stations, 
recorded  the  larger  surface  wave  periods.  This  factor  could  be 
related  to  the  instrumental  response  resulting  from  the  different 
galvanometers  used  in  the  seismograph  system.  Due  to  the  sparse 
horizontal  component  magnitude  data,  only  the  LPZ  component 
magnitude  was  graphically  plotted. 

The  overall  average  magnitude  (Mg)  was  U.09  which  was,  as 
expected,  lower  than  the  average  body  wave  magnitude  (Mg)  of  5.21 
(Figure  Hi).  Beyond  20  degrees,  the  surface  wave  magnitude  was 
U.53,  approximately  0.5  of  a  magnitude  unit  lower  than  the  computed 
body  wave  magni  .ude  (Mg)  for  distances  greater  than  16  degrees. 

The  characteristic  pronounced  shift  of  body  wave  magnitudes  with 
increasing  distance  was  not  reflected  by  the  surface  wave  magni¬ 
tudes.  The  magnitude  of  the  surface  waves  appeared  to  be  more 
evenly  distributed,  although  dispersed,  throughout  the  distance 
range.  As  noted  with  body  wave  magnitudes  beyond  16  degrees, 
there  was  a  corresponding  shift  of  surface  wave  magnitudes  in 
the  17  degrees  to  20  degrees  distance  range. 

The  surface  wave  magnitudes  (Figure  :.h)  do  not  appear  to 
correlate  with  azimuth,  wave  period,  or  region.  However,  they 
do  appear  to  be  distance  dependent,  indicating  a  review  of  the 
attenuation  function  may  be  necessary. 

Magnitudes  calculated  for  QASBUGGY  by  long-  and  short- 
period  instruments  at  the  same  station  generally  disagreed. 
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Figure  15  shows  the  Mg  -  Mg  (observed  at  the  same  station)  plotted 
versus  distance,  and  Figure  16  exhibits  similar  pairs  of  Mg  -  Mg 
measurements.  The  relatively  greater  range  of  values  shown  in 
Figure  l1?  nrobably  reflects  the  distance  dependence  of  both  Mg 
and  Mg.  The  obvious  improvement,  when  Mg  and  Mg  are  compared, 
indicates  a  fairly  consistent  relation  between  the  two  estimates 
of  magnitude.  In  general,  differences  between  the  surface  wave 
and  body  wave  magnitudes  decrease  with  increasing  distance. 

Figures  17  and  18  are  plots  of  Mg  vs  Mg  and  Mg  vs  Mg, 
respectively,  where  each  point  refers  to  a  pair  of  values  at 
the  same  station.  The  straight  line  segments,  which  were  devel- 
oped  by  Qutenberg  and  Richter  (1956),  and  by  Romney  (1963), 
represent  the  relationship  between  body  and  surface  wave  magnitudes 
of  earthquakes.  Again,  more  consistency  is  observed  when  Mg 
values  are  compared  with  Mg.  Perhaps  the  most  significant  in¬ 
formation  which  these  last  two  graphs  disclose  is  the  high  degree 
of  scatter  associated  with  pairs  of  magnitude  values  related 
to  the  same  source,  a  simple  compression,  as  compared  at  the 
same  station. 

A  general  comparison  concerning  Mg  versus  Mg  shows  that 
the  relative  excitation  of  long-period  surface  waxes  from 
explosions  is  less  than  that  from  shallow  earthquakes  since 
neither  WWNSS,  LRSM,  nor  Gutenberg  and  Richter  (1956)  curves 
appear  to  fit  the  GASHJGOY  data. 
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A  COMPARISON  BETWEEN  GASHJGGY  AND  AN  EARTHQUAKE 

On  January  23,  1966  at  Olh  56m  38s  U.T.,  an  earthquake  occurred 
at  36.96°N  and  106.90°W  which  placed  the  epicenter  circa  UO  kilo¬ 
meters  northeast  of  the  GASHJGGY  test  site.  As  both  epicenters  were 
essentially  in  the  same  place,  an  unusual  opportunity  was  provided 
for  a  comparison'  of  magnitudes,  both  Mg  and  Mg,  from  the  two  events. 

The  Qutenberg  and  Richter  (1956)  magnitude  form  la  was  used 
10  compute  the  body  wave  magnitudes,  whereas  the  surface  wave  magni¬ 
tudes  were  computed  by  the  formula  given  on  page  lb  in  this  report. 

Two  procedures  were  used  to  obtain  the  average  body  and  surface 
wave  magnitudes  for  the  two  events:  (1)  total  Mg  and  Mg  magnitude 
data}  (2)  ten  stations  which-  recorded  Mg  and  Mg  magnitudes  from 
GASHJGGY  and  from  the  earthquake.  The  total  data  resulted  in  an 
average  Mg  separation  of  0.11  unit  of  magnitude  between  the  two 
events  GASHJGGY:  5.21  (129  stations)}  Earthquake:  5.10  (3b 
stations  )J  .  The  average  body  wave  magnitude  for  the  ten-station 
network  differed  by  only  0.20  unit  (GASHJGGY  5.30}  Earthquake  5*10). 

The  average  surface  wave  magnitude  for  the  total  GASHJGGY  data  . 
was  U.09  (U2  stations)  and  for  the  earthquake  the  Mg  was  5*U2  (29 
stations).  This  is  a  difference  between  body  and  surface  wave 
magnitudes  of  1.12  units  for  GASHJGGY  and  for  the  earthquake  0.32  unit. 
A  comparison  of  the  average  surface  wave  magnitude  for  the  ten  stations 
shows  b.ll  for  GASHJGGY  and  5.37  for  the  earthquake.  Die  body  minus 


-  17  - 


surface  wave  magnitude  difference  for  GASBUGGY  was  1.19  magnitude 
units  as  compared  with  the  difference  for  the  earthquake  of  0.27 
unit. 

Figure  19  presents  a  plot  of  Mg  versus  Mg  for  the  earthquake 
where  each  point  refers  to  a  pair  of  values  at  the  same  station. 
Particularly  n9ticeable  in  Figure  19  is  the  apparent  increase 
or  decrease  of  Ms  ir  relation  to  Mg  as  compared  with  Figure  17 
at  similar  station; . 

Figure  20  illustrates  Mg  -  Mg  measurements,  which  were 
observed  at  the  same  station,  plotted  versus  distance  for  the 
earthquake.  The  majority  of  the  earthquake  data  are  positive 
which  is  in  contrast  to  the  negative  values  shown  by  the  GASBUGGY 
data  (Figure  15).  The  consistent  relationship  of  the  Mg  -  Mg 
values  is  noteworthy  for  similar  stations,  which  recorded  GAS&UGGY 
and  +v  earthquake,  in  the  h  degrees  to  16  degrees  distance 
range. 

The  difference  in  body  and  surface  wave  magnitude  is  one  method 
of  differentiating  between  a  shallow  earthquake  and  an  explosion. 

It  Jhas  been  noted  that  if  the  difference  is  approximately  one 
magnitude  unit  lower,  an  explosion  source  may  be  suspected.  The 
significant  difference  between  body  and  surface  wave  magnitudes, 
when  GASBUGGY  and  the  earthquake  are  compared,  tends  to  lend 
validity  to  the  earlier  observation. 
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HYPOC  ENTER  COMPUTATIONS 


The  results  of  the  computer  solutions  which  were  carried 
out  with  the  C&GS  hypocenter  program  (Ehgdahl  and  Ginst,  1966), 
are  shown  in  Table  The  computations  were  performed  with  the 
use  of  the  standard  travel-time  tables  of  Jeffreys-Bullen  (19U0) 
and  the  recently  published  Herrin  tables  (1968).  The  first  line 
in  Table  U  gives  the  known  coordinates  of  the  shot,  and  the 
remaining  lines  refer  to  computed  solutions. 

The  results  indicate  that  epicenters  can  be  located  within 
an  accuracy  of  a  few  kilometers  when  the  data  include  arrivals 
at  local  stations;  i.e.,  observations  of  the  "direct  wave".  The 
importance  of  having  such  nearby  observations  cannot  be  overstressed 
as  they  limit  the  error  in  location  and  are  essential  in  fixing 
the  depth  of  focus. 

As  has  been  noted  previously  (Carder,  et.al>,  1966),  the 
Jeffreys-Bullen  curve  is  a  good  first  approximation  of  regional 
travel  times  in  the  Western  mountain  environment.  An  examination 
of  residuals  from  all  computations  shows  that,  in  general,  they 
are  azimuthally  dependent,  and  that  they  indicate  particularly 
poor  agreement  with  the  standard  curves  in  the  distance  range 
of  10  to  20  degrees.  This  scatter  is  believed  due  to  the  obser¬ 
vation  that,  depending  on  signal -to -noise  conditions,  the  first 
identifiable  arrival  in  this  range  may  be  associated  with  the 
P7.9,  or  *io  5  Wsec  refractor. 


Nearly  all  observations  were  late  with  respect  to  the  1968 
curves,  indicating  velocities  beneath  the  GASEJJGGY  site  are 
relatively  slow  as  has  been  suggested  by  Herrinvand  Taggart's 
(1968)  seismic  delay  contours. 

HYPOC ENTER  RELOCATION 

An  attempt  to  refine  the  computed  hypoc enter  was  also  made 
by  adjusting  the  GASBUGGY  arrival  times  relative  to  the  residuals 
observed  from  GNOME  at  22  stations  which  recorded  both  events. 

The  results  were  spurious.  Using  a  restrained  depth  of  0  kilo¬ 
meters,  the  location  was  36.738°  N.,  107.1*1*7°  W. ,  and  a  location 
representing  unrestrained  parameters  was  36.720°  N.,  107.1*59°  W. 

The  shift  of  latitude  and  longitude  toward  the  northwest,  in 
kilometers,  from  the  known^coordinates  was  6.66  north,  21.1*1  west, 
and  for  the  latter  solution,  the  shift  was  l*. 66  north,  22.1*8  west. 
Apparently,  at  regional  distances  station  corrections  are  dependent 
on  azimuth  and  distance. 
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SUMMARY  OF  CONCLUSIONS 


1.  The  recorded  a  iplitudes  of  Pn  and  P  at  similar  distances 
were  equal,  or  the  amplitudes  were  slightly  higher  to  the  west 

of  GASBUGGY.  Amplitudes  of  Pn  and  P  appeared  to  be  higher 
where  the  first  arrival  was  early,  and  smaller  where  a  late 
arrival  was  recorded. 

2.  The  arrival  times  for  all  stations  appeared  to  follow 
the  Jeffreys-Bullen  curve  closely.  The  stations  to  the  east  of 
GASBUGGY  were  early,  but  not  as  early  as  found  in  GNOME.  Stations 
to  the  west  of  GASBUGGY  tended  to  be  early,  whereas  GNGKE  showed 
more  late  arrivals.  Beyond  20  degrees,  the  arrival  times  were 
about  1.0  sec  earlier  than  the  travel  time  predicted  by  the 
Jeffreys-Bullen  table. 

3.  Arrival  times  and  amplitudes  of  Pn  and  P  appeared  to  be 
correlated  with  physiographic  province. 

U.  Pg  and  Lg  recorded  amplitudes  were  higher  to  the  east 
than  to  the  west  of  GASBUGOY. 

5.  LQ  waves  were  recorded  to  the  east  and  to  the  west  of 
GASBUGGY.  The  wave  attenuation  rate  for  the  eastern  stations 
appeared  to  be  more  rapid  than  for  the  stations  to  the  west. 

6.  The  overall  average  magnitude,  when  computed  by  the 
method  given  by  Gutenberg  and  Richter  (195>6),  was  J>. 21 .  Using 
the  method  suggested  by  Evemden  (1967),  an  average  magnitude 
was  computed  to  be  U .73 • 


7.  The  recorded  data  confirm  the  Evernden  (196?)  apparent 


velocities,  which  indicate  the  absence  of  the  ^  velocity 
refractor  beneath  the  Great  Plains  to  the  east.  With  the  ex¬ 
ception  of  two  stations  to  the  east,  the  arrival  with  the  P, 
apparent  velocity  was  not  observed  east  of  the  Rocky  Mountain 
FVont. 


7.9 


8.  Rayleigh  waves  were  distinctly  recorded  at  large 
distances  from  GASBUGGY,  whereas  GNOME  recorded  LR  waves  out  to 
a  distance  of  500  kilometers.  IVenty  second  surface  waves  were 
not  recorded  from  GASBUGGY.  The  surface  wave  magnitude  (Mg)  was 
U.09  «3  compared  with  a  Mb  and  Mg  of  5-21  and  li . 73 ,  respectively. 
The  observation  that  Ms  is  wore  than  one  magnitude  unit  less  than 
Mb  agrees  with  previously  published  results  in  discerning  earth¬ 
quakes  from  explosions. 

9.  The  attempt  to  refine  the  GASBUGGY  location,  by  cor¬ 
recting  the  arrival  times  relative  to  the  observed  residuals  from 
GNOME,  indicates  that  location  accuracy  is  dependent  on  recording 
station  distance  and  azimuthal  control. 


« 
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TABLE  1#  ••  GAS  BUGGY  Event ,  Arrival  Tines  end  Amplitudes  of  Principal  Phases 
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Obs.  Arrival  Tlae  Period  Maxima 

Distance  Con-  1900  Hrs.  T  A^>lltade  Magnitude  J.B. 
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Obs.  Arrival  Tim  Period  Maxima 

Distance  Com-  1900  Hr*.  T  A^lltude  Magnitude  J.B. 

CODE _ Station _ Degrees  ta»  ponent  Phase  Min. _ See. _ Sec.  A/T  (an/ sec)  Mp  Mr  Residual 
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Obs.  Arrival  Tine  Period  Maxlsum 

Distance  Com-  1900  Hrs.  T  Altitude  Magnitude  J.B. 

CODE  Station  Degrees  km  ponent  Riase  Min.  •  Sec.  Sec.  A/T  (m/i/sec)  Mp  Mp  Residual  Remarks 

WH2T1C  White  Horse  29. 79  3311*  SPZ  eP  36  (09.7)  0.7  U.l  U.2  -0.8 
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FIG.  2  DISTANT  SEISMOGRAPHIC  STATION',  RECORDING  SIGNALS  FROM  GASBUGGY. 
THE  SMOOTH  CURVE  IS  AT  A  DISTANCE  OF  100°  FROM  GASBUGGY  SITE. 
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Appendix  2  --  Q  Tables  for  Specific  Apparent  Velocities 
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